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Abstract. We have performed a set of 1 1 three-dimensional magnetohydrodynamical core-collapse supernova 
simulations in order to investigate the dependencies of the gravitational wave signal on the progenitor's initial 
conditions. We study the effects of the initial central angular velocity and different variants of neutrino transport. 
Our models are started up from a 15M G progenitor and incorporate an effective general relativistic gravitational 
potential and a finite temperature nuclear equation of state. Furthermore, the electron flavour neutrino transport 
is tracked by efficient algorithms for the radiative transfer of massless fermions. We find that non- and slowly 
rotating models show gravitational wave emission due to prompt- and lepton driven convection that reveals 
details about the hydrodynamical state of the fluid inside the protoneutron stars. Furthermore we show that 
protoneutron stars can become dynamically unstable to rotational instabilities at T/\W\ values as low as ~ 2% at 
core bounce. We point out that the inclusion of deleptonization during the postbounce phase is very important 
for the quantitative GW prediction, as it enhances the absolute values of the gravitational wave trains up to a 
factor of ten with respect to a lepton-conserving treatment. 



PACS numbers: 04.30.Db, 95.30.Qd, 97.6().Bw 



1. Introduction 



Stars in the mass range 8M < M < 4OM end their lives in a core-collapse supernova (CCSN). However, 
at present the fundamental explosion mechanism, which causes a star to lose its envelope by a yet uncertain 
combination of factors including neutrino heating, rotation, hydrodynamical instabilities, core g-mode 
oscillations and magnetic fields, is still under debate (for a review, see e.g. dl lb). As strong indications both 
from theory and observations exist that CCSNe show aspherical, multidimensional features there 
is a reasonable hope that a small amount of the released binding energy will also be emitted as gravitational 
waves (GWs), thus delivering us first-hand information about the dynamics and the state of matter at the 
centre of the star. GW emission from CCSNe were suggested to arise from i) axisymmetric rotational 
core collapse and bounce ( 51 ; "nt 5 ; HI) ii) p rompt-, neutrino-driven postbounce convection and anisotropic 
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) . iii) protoneutron star (PNS) g-mode oscillations d35l) and 
iv) nonaxisymmetric rotational instabilities (Hit 48; 39; 40|; 38; 37; 42; 43) For recent reviews with a more 
complete list of references, see (2 34 )■ However, only i) can be considered as being well understood as far 
as the physics of the collapse is concerned, since only theses models incorporate all relevant input physics 
known at present ®) (there are, though, still large uncertainties with respect to the progenitor star, e.g. 
rotation profiles, magnetic fields, and inhomogenities from convection). The prediction of all other suggested 
emission scenarios (ii-iv) still neglect, to a certain extent, dominant physics features due to the diversity and 
complexity of the CCSN problem on the one hand side and restrictions of available computer power on 
the other side. Hence, the computational resources were so far either spent on highly accurate neutrino 
transport (e.g. (21; 36[ 2^)) while neglecting other physical degrees of freedom such as magnetic fields, or 
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focus on a general relativistic treatment and/or 3D fluid effects such as accretion funnels, rotation rate and 
convection, but approximate or even neglect the important micro physics. Only recently have detailed 3D 
computer models of CCSN become feasible with the emerging power of tens of thousands CPUs unified in 
a single supercomputer. Such detailed simulations are absolutely indispensable for the following reasons: 
a) GW astronomy requires not only very sensitive detectors, but also depends on extensive data processing 
of the detector output on the basis of reliable GW estimates (Q]). b) The temperatures and densities inside 
a supernova core exceed the range that is easily accessible by terrestrial experiments. Thus, it will be 
impossible for the foreseeable future to construct a unique finite temperature equation of state (EoS) for hot 
and dense matter based on experimentally verified data. Therefore, models with different parameter settings 
must be run and their computed wave form output then can be compared with actual detector data. Hence, 
modelling will bridge the gap between theory and measurement and allowing the use of use of CCSNe 
as laboratory for exotic nuclear and particle physics d25l) . In this paper, we will present the GW analysis 
of a set of 1 1 three-dimensional ideal magnetohydrodynamical (MHD) core-collapse simulations. We will 
focus our study on the imprint of 3D nonaxisymmetric features onto the GW signature. Our calculations 
include presupernova models from stellar evolution calculations, a finite-temperature nuclear EoS and 
a computationally efficient treatment of deleptonization during the collapse phase. General relativistic 
corrections to the spherically symmetric Newtonian gravitational potential are taken into account. Moreover, 
while several models incorporate long-term neutrino physics by means of a leakage scheme, we also present 
the first results of a model which includes a neutrino transport approximation in the postbounce phase that 
takes into account both neutrino heating and cooling. As for the progenitor, we systematically investigate 
the effects of the spatial grid resolution, the neutrino transport physics and the precollapse rotation rate with 
respect to its influence on the nonaxisymmetric matter dynamics. 

This paper is organised as follows. In section 2 we briefly describe the initial model configurations and 
the numerical techniques employed for their temporal evolution. Furthermore, we review the tools used for 
the GW and data analysis. Section 3 collects the results of our simulations. Finally, section 4 contains our 
conclusions and an outlook of our future research. 



1.1. Description of the magnetohydrodynamical models 



For the 3D Newtonian ideal MHD CCSN simulations presented in this paper, we use the FISH code (Il2l) . 
The gravitational potential is calculated via a spherically symmetric mass integration that includes radial 
general relativistic corrections (0). The 3D computational domain consists of a central cube of either 600 
or 1000 3 cells, treated in equidistant Cartesian coordinates with a grid spacing of 1km or 0.6km. It is, as 
explained in d42l) . embedded in a larger spherically symmetric computational domain that is treated by the 
time-implicit hydrodynamics code 'Agile' (1231) . Closure for the MHD equations is obtained by the softest 
version of the finite-temperature nuclear EoS of (1191) . The inclusion of neutrino physics is an essential 
ingredient of CCSNe simulations, as ~99% of the released binding energy is converted into neutrinos of all 
flavours. Their complex interactions with matter (e.g. dill) ) are believed to drive the supernova explosion 
dynamics in the outer layers as well as deleponizing the PNS to its compact final stage as a neutron star. As 
the Boltzmann neutrino transport equation can only be numerically solved in a complete form in spherical 
symmetry on today's supercomputers (30). our 3D simulations must rely on several feasible approximations 
which capture the dominant features of the neutrino physics. As for the treatment of the deleptonization 
during the collapse phase, we apply a simple and computationally efficient Y e vs. p parametrization scheme 
which is based on data from detailed ID radiation-hydrodynamics calculations ( 1240 . For this we use the 
results obtained with the Agile-Boltztrancode d22l) . including the above-mentioned EoS and the electron 
capture rates from Around core bounce, this scheme breaks down as it cannot model the neutronization 
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Table 1. Summary of the models' initial conditions and GW related quantities. fi c , [rads ] is the precollapse 
central angular velocity, while = T/\W\ is the ratio of rotational to gravitational energy. p c j, [10 14 gcm -3 ] is 
the maximum central density at the time of core bounce. Eaw [10 -9 M Q c 2 ] is the energy emitted as GWs. 
[Hz] denotes the peak frequency of the GW burst at bounce, while fyw [Hz] stands for the spectral peak from 
the narrow band emission caused by a low T/\W\ instability, tf [ms] is the time after core bounce when the 
simulation was stopped. 
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burst. After core bounce, the neutrino transport thus is tracked for several models via a partial (i.e. a leakage 
scheme) or full implementation of the isotropic diffusion source approximation scheme (IDSA, J^t)). The 
IDSA decomposes the distribution function / of neutrinos into two components, a trapped component f 
and a streaming component f s , representing neutrinos of a given species and energy which find the local 
zone opaque or transparent, respectively. The total distribution function is the sum of the two components, 
f - f + f ■ The two components are evolved using separate numerical techniques, coupled by a diffusion 
source term E. The source term £ converts trapped into streaming particles and vice versa. We determine 
it from the requirement that the temporal change of /' has to reproduce the diffusion limit in the limit of 
small mean free path. Note that our leakage scheme significantly overestimates the deleptonization in and 
around the neutrino sphere region, as it neglects any absorption of transported neutrinos by discarding the 
streaming component (/ s = 0). The presupernova stellar models stem from Newtonian ID stellar evolution 
calculations and hence may not cover all possible states prior to the collapse of a multidimensional star. 
Therefore we construct the initial conditions of our simulations by a parametric approach. We employ a 
solar-metallicity 15M progenitor of KEdti . and set it into rotation according to a shell-type rotation law of 
(0) with a shellular quadratic cutoff at 500km radius. The initial magnitude of the magnetic field strength for 
all models is fixed at values suggested by 43). 



1.2. Gravitational Wave extraction 

We employ the Newtonian quadrupole formula in the first-moment of momentum density formulation (HI) 
to extract the GWs from our simulation data. Note that the quadrupole formula is not gauge invariant and 
only valid in the Newtonian slow-motion limit ( 31 ). Nevertheless, it was shown by d44h in comparative tests 
to work sufficiently well compared to more sophisticated methods, as it preserves phase while being off in 
amplitude by ~10%. 
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Figure 1. Left: Time evolution of the GW + polarization for a spectator located at the polar axis (Models Rl#s, 
RO/dsa and R1l) Right: Corresponding spectral energy distribution of models RO/osa and Rl/, at a distance of 
lOkpc compared with the LIGO strain sensitivity and the planned performance of Advanced LIGO (45 ). Optimal 
orientation between source and detector is assumed. 




x [ie7 cm] Mass [Ms] 



Figure 2. Left: Model RQidsa's specific entropy distribution [kg/baryon] (left side) and electron fraction Y e 
(right side) 50ms after core bounce. The innermost 600 2 km 2 in the x-y plane are displayed. The entropy color 
bar scales from (blue) to 12 (red). The Y e color bar accounts for values from (red) to 0.5 (blue) Right: 
Comparison of the spherically averaged Y e profiles of models ROl (circled line, 'LEAK'), RO/dsa (solid line, 
'IDS A') with the spherically symmetric model G15 (dashed line, 'Boltzmann') as a function of the enclosed 
mass at 5ms after bounce. Model G15 is based on general relativistic three-flavour neutrino transport (52). 



2. Results 

2.1. Non- or slowly rotating core collapse 

In order to study the influence of neutrino transport on the stochastic matter dynamics in the early supernova 
stages (t < 100ms after bounce), without having other different physical parameters interfering, we 
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carried out three simulations: RO (purely hydrodynamical postbounce evolution), Rl^ (includes a leakage 
scheme) and ROidsa (incorporates both neutrino cooling and heating). Non- and slowly rotating progenitors 
(Q c , < 0.3rads~' in our model set) all undergo quasi-spherically symmetric core collapse. As the emission 
of GWs intrinsically depends on dynamical processes that deviate from spherical symmetry, the collapse 
phase therefore does not provide any kind of signal, as shown in Fig |l(a)| for t—tb < 0. However, subsequent 
pressure-dominated core bounce, where the collapse is halted due to the stiffening of the EoS at nuclear 
density p mlc » 2 x 10 14 gcirT 3 , launches a shock wave that plows through the infalling layers, leaving behind 
a negative entropy gradient. Moreover, as soon as the shock breaks through the neutrino sphere ~ 5ms 
after bounce, the immediate burst of electron neutrinos causes a negative lepton gradient at the edge of 
the PNS. The combination of these two gradients form a convectively unstable region according to the 



Schwarzschild-Ledoux criterion d 1 8c 1491) . which in turn induces a GW burst due to this so-called 'prompt' 
convection. A detailed comparison of the models RO, RIl, RO/as^ shows that all of them follow a similar 
dynamical behaviour until about 20ms after bounce. At this stage, aspherities leading to GW emission are 
predominantly driven by the negative entropy gradient and not by the lepton gradient. Hence, the wave trains 
of all three models, which are based on stochastic processes, fit each other relatively well both in amplitude 
(several cm) and spectra (~ 150 - 500Hz). However, 'prompt' convection depends, as it was pointed out 



by d34l) . not on the negative entropy gradient alone, but also on numerical seed perturbations which are 
introduced by the choice of the computational grid. Hence, in order to test the dependence of our findings 
on the spatial resolution, we carried out model Rl//«. This better resolved simulation shows considerably 
smaller seed perturbations around t - t j ~ 0, as grid alignment effects are better suppressed at core bounce; 
hence prompt convection then is much weaker and a smaller GW amplitude (~ 50%) is emitted, as shown in 
Fig |l(a)| However, better numerical resolution also leads to less numerical dissipation in the system, which 
eases the dynamical effects that follow. Thus, we find for ~ 10ms < t < 20ms considerably stronger GW 
emission in RIhr compared to the 1km resolved models, as indicated in Fig |l(a)| The three representative 
simulation results diverge strongly in the later postbounce evolution (t > 20ms). Convective overturn causes 
a smoothing of the entropy gradient. As a result, the GW amplitude in the hydrodynamical model R0 
quickly decays (t < 30ms after bounce) and is not revived during the later evolution. On the other hand, the 
negative radial lepton gradient (see Figs |2(a)| and |2(b)| which is caused by the neutronization burst and the 
subsequent deleptonization, which we model only in R1/, and RO/asa, now starts to drive convection inside 
the PNS. For the latter models, the so-called PNS convection (0) exhibits similar maximum amplitudes of 
~l-2cm (Fig |l(a)[ ), while differing from each other strongly in the corresponding spectra, as displayed in 
Fig |l(b)| RIl's spectrum peaks between ~600 -1000Hz, while RO/asa's frequency band peaks at values as 
low as -100Hz. This affects the total energy E cw emitted (O(lO" lo )M c 2 vs. 6>(lO" n )M c 2 , see TabfJ}, 
being one order of magnitude higher for RIl due to dEcw/df °c f 2 . We found the key controlling factor 
of this behaviour to be the radial location of the convectively unstable zones and the related dynamical 
characteristical timescales f ( / v „ involved. If we use as rough estimate t^ v „ ~ A r /c7, 0, and apply typical 
values for the models ROidsa and RIl, we confirm the obtained values. Furthermore, our leakage scheme 
significantly overestimates neutrino cooling processes, as one can see in Fig |2(b)| There, the convectivly 
unstable layer is extended to radii above nuclear densities, where matter still is opaque for neutrinos and 
where the local speed of sound assumes values far larger than in the case of model RO/asa- Hence, the 
dynamical timescale of RIl is considerably shorter and the spectral distribution is peaked at higher values. 
When comparing the results found for model ROidsa with a very recent 2D study of (f29h . where they carried 
out one simulation (cf. their model M15LS-2D) with comparable input physics (same 15M Q progenitor; 
same underlying finite-temperature EoS) and a very sophisticated neutrino transport scheme, we find very 
good agreement both in the amplitudes and frequencies. Hence we conclude that the primary ingredient 

$ cj = 1 /A, J r c s (f)dr is the radially averaged sound speed of a convectively unstable layer with a radial extension of A r . 



Gravitational waves from supernova matter 



6 




Figure 3. Left: Vorticity's z-component w, = (V X v) z of model R5 L in the equatorial plane, 29ms after bounce, 
showing a dominant m=l mode. The innermost 300 2 km 2 are displayed, and the color is encoded in units of 
[s -1 ], ranging from -5000 (white) to 5000 (black) Right: GW polarizations + and X for models R5^ and R5 as 
seen from an observer along the polar axis. Strong nonaxisymmetric dynamics with m=2 components develop 
right after core bounce. Hence, the two polarizations are shifted by a quarter cycle, as one could expect from 
GWs emitted by a spinning bar. 



for supernova simulations which attempt a quantitative prediction of GWs from 'prompt' and early PNS 
convection (t < 100ms after bounce) is the accurate radial location and size of convectively unstable layers. 
It defines the dynamical behaviour and timescale of overturning matter in this early supernova stage. 

2.2. Rotational core collapse & nonaxisymmetric instability at low T/\W\ 

Recently it has been argued based on numerical simulations of equilibrium neutron star models or full 
core-collapse simulations that differentially rotating PNS can be subject to non-axiymmetric rotational 
instabilities (see Fig |3(a)[ ) at /3 values (=T/\W\, the ratio of rotational to gravitational energy) far below 
the ones known from the classical dynamical bar mode instability with a threshold of = 27%, or the 
secular instability, which is triggered at p sec ~ 14% (0), leading to strong, narrow-band GW emission, as 



displayed in Fig |3(b)| (|3t l37t 1421) . At present little is known about the true nature of the so-called low T/\W\ 
instability. Previous work has so far failed to establish (for example) an analytical instability criterion, as was 
pointed out by ( I341) . We addressed two relevant questions regarding the so-called 'low T/\W\' instability in 
the context of stellar core collapse: i) Which is the minimum /3 value required in self-consistent core-collapse 
simulation to trigger the onset of the instability? This is important to know since most stars which undergo 
a core collapse rotate only slowly (0); furthermore, it was pointed out by (0) that even fast rotating PNS can 
never accrete enough angular momentum to reach the /3d y „ value required for the onset of the classical bar 
mode instability, ii) How does the inclusion of deleptonization in the postbounce phase quantitatively alter 
the GW signal? So far, 3D models have not included spectral neutrino physics in the postbounce phase. To 
study i), we systematically change the rotation rate while keeping the other model parameters fixed. The 
minimum T/\W\ value we found in our parameter range to trigger the instability was ~ 2.3% at core 
bounce (model R3), which is considerably lower than seen in previous studies ( d37l) found fa, ~ 9%, while 



d42l) found fa, ~ 5%). Furthermore, we find that centrifugal forces set a limit to the maximum frequency of 
the GW signal around ~ 900Hz. The faster the initial rotation rate, the stronger the influence of centrifugal 
forces, which slow down the postbounce advection of angular momentum onto the PNS. The result is a 
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slower rotation rate, a lower pattern speed and thus GW emission at lower frequencies (see TabQ}. In order 
to address ii), we carried out 'leakage' model R5l. This model shows 5-10 x larger maximum amplitudes 
due to the nonaxisymmetric dynamics compared to its hydrodynamical counterpart R5 that neglects neutrino 
cooling (see Fig. |3(b)| i. This suggests that the treatment of postbounce neutrino cooling plays an important 
role when it comes to the quantitative prediction of GW signals from a low /3 instability. The neutrino cooling 
during the postbounce phase leads to a more compact PNS with a shorter dynamical timescale compared to 
the purely hydrodynamical treatment. This in turn is reflected in the dynamical evolution. The shock wave 
stalls at considerably smaller radii and becomes more quickly unstable to azimuthal fluid modes. Since 
there is much more matter in the unstable region of this model, the unstable modes grow faster, causing the 
emission of much more powerful GWs. However, we again point out that our leakage scheme overestimates 
the compactification of the PNS due to neutrino cooling. The 'reality' for the strength of GW emission 
therefore should lay in between the results from the pure hydrodynamical- and leakage treatment. 



3. Summary and outlook 

We have presented the GW signature of eleven 3D core-collapse simulations with respect to variations 
in the spatial grid resolution, the underlying neutrino transport physics and the initial rotation rate. Our 
results show that in case of non- and slowly rotating models the GWs emitted during the first 20ms 
after bounce are predominantly due to entropy driven 'prompt' convection. It turns out that the crucial 
parameter to study this stochastic phenomenon is the choice of the spatial resolution and not the inclusion 
of a neutrino transport scheme. This parameter has a twofold effect: Firstly, it governs the influence of 
numerical noise, since a better resolution leads to lower numerical seed perturbations and thus smaller 
grid alignment effects. Therefore, the GW amplitude right at core bounce is smaller for higher spatial 
resolution. Secondly, it enhances the ability to follow dynamical features, as better numerical resolution 
causes less numerical dissipation in the system, which eases the dynamical effects which follow, leading 
to larger GW amplitudes after the core bounce compared to less resolved models. The lepton driven 
convection is the central engine for the later dynamical postbounce evolution of the PNS (t > 20ms) and 
hence the GW emission. Our findings and comparisons with state of the art 2D simulations of J29I) suggest 
that the radial location and size of the convectively unstable layers are the key controlling factor for the 
outcome of the GW prediction, as they define the timescale and the dynamical behaviour of the overturning 
matter. Here we find a large sensitivity to the numerical approach of the neutrino transport scheme. In our 
rotational core-collapse simulations, nonaxisymmetric dynamics develops for models with a rotation rate of 
Pb <: 2.3% at core bounce. Beyond this value, which is considerably lower than found in previous studies 



(e.g. (I37t 1421) 1. all models become subject to the 'low T/\W\' instability of dominant m=l or m=2 character 
within several ms after bounce (Fig |3(a)[ ). The fact that the effectively measured GW amplitude scales with 
the number of GW cycles N as h e ff oc h -y^V (47) suggests that the detection of such a signal is tremendously 
enhanced. Moreover, we point out that the inclusion of deleptonization during the postbounce phase causes 
a compactification of the PNS which enhances the absolute values of the GW amplitudes up to a factor of 
ten with respect to a lepton-conserving treatment. 

The major limitation of our code now is in the monopole treatment of gravity, since it cannot account 
for spiral structures, which could be reflected in GW. We are currently working on the improvement of 
this issue. Furthermore, the IDSA includes at present only the dominant reactions relevant to the neutrino 
transport problem (see J26I) for details). Future upgrades will also include contributions from electron- 
neutrino scattering, which are indispensable during the collapse phase. The inclusion of this reaction will 
also make the cumbersome switch of the parametrization to the IDSA at bounce obsolete. Finally, we work 
on the inclusion of p. and r neutrinos, which are very important for the cooling of the PNS to its final stage 
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